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[57] ABSTRACT 

The instant invention provides a method and apparatus for 
use in laser induced dynamic light scattering which attenu- 
ates the multiple scattering component in favor of the single 
scattering component. The preferred apparatus utilizes two 
light detectors that are spatially and/or angularly separated 
and which simultaneously record the speckle pattern from a 
single sample. The recorded patterns from the two detectors 
are then cross correlated in time to produce one point on a 
composite single/multiple scattering function curve. By col- 
lecting and analyzing cross correlation measurements that 
have been taken at a plurality of different spatial/angular 
positions, the signal representative of single scattering may 
be differentiated from the signal representative of multiple 
scattering, and a near optimum detector separation angle for 
use in taking future measurements may be determined. 

18 Claims, 4 Drawing Sheets 
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METHOD AND APPARATUS FOR FIBER 
OPTIC MULTIPLE SCATTERING 
SUPPRESSION 

RELATED APPLICATIONS 

This application claims the benefit of U.S. Provisional 
Application No. 60/059,619, which application was filed 
with the Patent and Trademark Office on Sep. 23, 1997, the 
disclosure of which is incorporated herein by reference. 

STATEMENT REGARDING FEDERALLY 
SPONSORED RESEARCH 

This invention was made with Government support under 
contract Number NAG3-1624 awarded by NASA and Grant 
Number DMR-9501865 awarded by NSF. The Government 
has certain rights in this invention. 

BACKGROUND OF THE INVENTION 

FIELD OF THE INVENTION 

This invention relates generally to determining the physi- 
cal properties of materials through the use of dynamic light 
scattering, and, more specifically, to a technique for detect- 
ing and evaluating singly scattered speckle in concentrated 
solutions by suppressing the multiply scattered contribution. 

BACKGROUND 

When light impinges on matter, a portion of that light is 
scattered or reflected away. The frequency shifts, the angular 
distribution, the polarization, and the intensity of the scat- 
tered light are determined by the size, shape, and molecular 
interactions of the particles in the scattering material. Thus, 
from the light-scattering characteristics of a given system it 
is possible, with the aid of electrodynamics and the theory 
of time -dependent statistical mechanics, to obtain informa- 
tion about the structure and molecular dynamics of the 
scattering medium. 

Dynamic light scattering exploits the time dependent 
coherence loss of scattered light to explore the movement of 
the scatterers. Its applications range from measurements of 
molecules or particles smaller than the wavelength of the 
probing light [1] to the much larger structures of, for 
example, foam [2]. Small scatterers are easily analyzed 
through the single scattering of photons. Single scattering 
occurs when each molecule or particle is exposed to essen- 
tially the same incident light as occurs in relatively dilute 
solutions. 

Diffusing Wave Spectroscopy [3] and Correlation Trans- 
fer Theory [4] make the investigation of more concentrated 
colloidal samples possible. Both theories, however, presume 
a knowledge of the single scattering correlation function. 
This function is often unknown in practical multiple scat- 
tering applications, where the incident light may not equally 
expose all particles or where certain particles may be 
exposed to light scattered by other particles. Further, even 
though an equivalent single scattering decay time can be 
obtained by using the transmitted intensity relative to a 
known optical thickness reference, where the standard func- 
tional form for the single scattering correlation function as 
an exponential is known, the form of the mean square 
displacement must still be assumed. It would therefore be 
very useful to measure the single scattering contribution 
inside a multiple scattering sample. 

Heretofore, two different schemes have been available for 
extracting some single scattering information from a mul- 
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tiple scattering sample. One such scheme is a fiber semi- 
backscattering technique that utilizes scattered photons from 
a small overlap region between an input and a detecting fiber 
[5]. In this approach, most of the scattering paths within the 
5 sample volume involve only single collisions with particles. 
However, this method does not yield a systematic way of 
determining the degree of single scattering for different 
samples. Additionally, in some cases the solution, or the 
particles present therein, can coat the detecting fiber 
10 (because the fiber may need to be inserted into the sample 
during measurement), thereby corrupting the received sig- 
nal. 

The other approach focuses on multiple scattering elimi- 
nation by cross correlating the signal of two different inci- 
15 dent wave lengths at the same scattered wave vector for each 
[6, 7] and is based on the multiple scattering elimination 
methods of Phillies [8, 9, 10]. Two -Color-Cross- Correlation 
detects only true single scattering. A drawback of this 
technique is that it requires extensive careful alignment of 
20 two input laser beams and two detectors. Additionally, the 
turbidity of the sample must be small enough to still exhibit 
significant single scattering. Thus, the technique is limited to 
small scattering concentrations that can be increased slightly 
by reducing the sample dimensions. 

25 It is accordingly an object of the present invention to 
analyze the single scattering of light in turbid media. 

It is a related object to provide an efficient means of 
suppressing the multiple scattering contribution from an 
illuminated turbid media so that true single scattering data 
may be obtained and evaluated. 

It is another object that the invention be cost effective 
from an equipment requirement standpoint, utilizing a single 
laser input beam and readily available component parts. 

35 Before proceeding to a description of the instant 
invention, however, it should be noted and remembered that 
the description of the invention which follows, together with 
the accompanying drawings, should not be construed as 
limiting the invention to the examples (or preferred 
40 embodiments) shown and described. This is so because 
those skilled in the art to which the invention pertains will 
be able to devise other forms of this invention within the 
ambit of the appended claims. 

SUMMARY OF THE INVENTION 

45 

Broadly speaking, the instant invention provides a method 
and apparatus for estimating single scattering functions - 
particularly in concentrated solutions. The instant method 
utilizes two light detectors that are spatially and/or angularly 
50 separated and which simultaneously record the speckle 
pattern from a single sample. The recorded patterns from the 
two detectors are then cross correlated at zero-lag to produce 
one point on a composite single/multiple scattering function 
(S/N) curve. By collecting and analyzing cross correlation 
55 measurements that have been taken at a plurality of different 
spatial/angular positions, the signal representative of single 
scattering may be differentiated from the signal representa- 
tive of multiple scattering, and a near optimum detector 
angle for use in taking future measurements may be deter- 
60 mined. 

By way of general background, when a laser beam is 
directed into a scattering medium and reflects from a surface 
a speckle pattern is often observed. This complex pattern 
results from interference of electromagnetic radiation which 
65 originated from a coherent source but which has followed 
different paths in reflecting or scattering to the detector. At 
some points, the total field reflecting from the surface will 
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add constructively and be bright, while at other points the 
total field will add destructively and be dark. If the laser 
beam is focused to a small region on the reflecting surface, 
the speckle size typically increases in dimension at the 
detector. This diffraction effect is analogous to that observed 5 
in single slit diffraction, where the diffraction pattern width 
increases as the slit width decreases. Further, the single 
scattered light, arising as it does within the laser beam, will 
tend to originate from a smaller region than the multiple 
scattered light which tends to be diffused throughout the 
sample. For this reason, at the detector the single scattered 
light will have a broad speckle field compared to the 
multiple scattered speckle. 

Turning now to a discussion of various aspects of the 
instant invention, according to one preferred embodiment 
there is provided a method, whereby it is possible to identify 15 
and separate the contribution of the multiple scattered light 
from the single scattered light based on the properties of the 
respective speckles. Since the single scattering speckle 
(which arises from inside of the incident laser beam) is 
correlated over a wider angular or spatial range than is the 20 
multiple scattering speckle (which might potentially origi- 
nate from anywhere within the sample), by using two 
detectors with an appropriately selected spatial (or angular) 
separation, it is possible to detect the single scattered speckle 
and attenuate or exclude the multiple scattered speckle 25 
contribution. At each angular setting, the light intensities 
received within the two detectors are cross correlated in time 
to produce a single numerical value for that detector angle. 

By collecting cross correlations at a number of different 
detector angles, a profile is created that can be used to 30 
determine the angular ranges over which multiple and single 
scattering predominate. From the individual cross correla- 
tions and the ensemble of angle and cross correlation pairs, 
it is additionally possible to estimate the optimum (or near 
optimum) detector angle to use in collecting future 35 
observations, as well as a variety of other useful quantities 
that pertain to the particles of the laser-illuminated sample. 

A preferred apparatus for use with the instant invention 
suppresses multiple light scattering by using two slightly- 
tilted detectors that are directed toward the same sample 40 
volume which has been illuminated by a narrow laser beam. 
Further, a mechanism is provided whereby the angle 
between the two detectors may be continuously adjusted so 
as to give readings at any selected separating angle. The 
apparatus further preferably utilizes single-mode fiber optic 45 
fibers as detectors, single mode fibers being almost ideal for 
use in detecting dynamic light scattering in dilute suspen- 
sions. Additionally, optical fiber GRIN (i.e., graded index) 
lenses are preferably used on the ends of the monomode 
fibers to improve the performance of the invention. 50 

As another preferred embodiment of the instant invention, 
there is provided an apparatus substantially as described 
above, but wherein a polarizing filter is placed between the 
sample and the detectors. In the preferred embodiment, the 
polarizing filter is placed in one of two orientations with 55 
respect to the incident laser beam’s polarization plane: it is 
either oriented to have the same polarization (parallel) or 
oriented to have a perpendicular polarization 
(perpendicular). Although both polarization components 
show multiple scattering contributions, only the parallel 60 
component contains evidence of the single scattering signal. 
Thus, the orientation of the polarizer may be used as part of 
a preferred embodiment to increase the sign al-to -noise ratio 
and eliminate some of the multiple scattering in the received 
signal. 65 

According to still another preferred embodiment, the 
collected ensemble of angle/cross correlation pairs can be 
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used with numerical curve fitting routines to provide an 
estimate of, by way of example, a two-cumulant represen- 
tation of the time decay of the intensity correlation function 
and, thus, an estimate of the particle size. When applied to 
the polarization embodiment of the instant invention, the 
curve fitting results for the two polarization components 
support the theory which calls for larger single scattering 
speckles in the parallel component, and smaller multiple 
scattering speckles in both the parallel and perpendicular 
components. 

Finally, in another preferred embodiment the methods 
described previously are extended to much higher concen- 
trations by detecting the scattering from particles close to the 
sample boundaries. When the instant apparatus detects the 
particles close to the boundaries of the sample container, 
there is no limit on the applicable optical thickness except 
for the exponential decay of intensity. There are always 
some multiple scattering contributions present that lead to 
some non-exponentiality. The main decay, though, is sur- 
prisingly accurate and can be used to measure the particle 
sizes of concentrated samples with ease. 

Abetter understanding of the present invention, its several 
aspects, and its objects and advantages will become apparent 
to those skilled in the art from the following detailed 
description, taken in conjunction with the attached drawings, 
wherein there is shown and described the preferred embodi- 
ment of the invention, simply by way of illustration of the 
best mode contemplated for carrying out the invention. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 is a diagram of an experimental setup for the 
preferred fiber optic multiple scattering suppression tech- 
nique. 

FIG. 2 is a diagram of the preferred detecting assembly. 

FIG. 3 is a representative plot of intercept versus angular 
tilt for the cross correlation of a multiple scattering sample. 

FIG. 4 is a representative plot of intercept versus angular 
tilt for the two polarizations of a multiple scattering sample. 

FIG. 5 is a representative plot of particle size (nm) versus 
angular tilt for the two polarizations of a non-single scat- 
tering sample. 

FIG. 6 is a representative plot of normalized second 
cumulant versus angular tilt for two polarizations of a 
non-single scattering sample. 

DETAILED DESCRIPTION OF THE 
PREFERRED EMBODIMENT 

Before explaining the present invention in detail, it is 
important to understand that the invention is not limited in 
its application to the details of the construction illustrated 
and the steps described herein. The invention is capable of 
other embodiments and of being practiced or carried out in 
a variety of ways. It is to be understood that the phraseology 
and terminology employed herein are for the purpose of 
description and not of limitation. 

THEORETICAL BACKGROUND 

Broadly speaking, the instant invention utilizes two detec- 
tors that are aimed at slightly different angles at a single 
laser-illuminated sample to record scattered light emanating 
therefrom. Measurements taken from the detectors over a 
range of detector angles are used to determine a best angular 
separation to use so as to observe predominantly the single 
scattered light component. 

By way of explanation, for a Gaussian random process, 
the intensity correlation function between the two detectors 
takes the following form [15]: 
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<i(7 lf Qi(7 2 ^ 2 )>=<i(7 1 ,r 1 )><i(7 2 ^ 2 )>[i+Iy(7 l9 T 2 ,h-t 2 )\ 2 l (i) 

In this equation, the function y( # ) is the second order 
complex degree of coherence, and I( r , t) is the light 

intensity at the point r at time t. The bracket notation (i.e., 
“<*>”) indicates that a time average or expected value is to 
be taken of the argument (or an ensemble average for a 
stationary process). In typical dynamic light scattering 

applications, there is a single detector ( r 1 = r 2 ) and the 
time dependence of the (auto) correlation function carries 
the desired information. However, of particular interest for 
purposes of the instant invention is the signal-to -noise ratio 
or the complex degree of coherence when t 1 =t 2 , for a similar 
signal measured by two detectors oriented at slightly differ- 
ent angles. If intensities originating at different regions in the 
scattering volume are spatially uncorrelated, then the van 
Cittert-Zernike theorem in the far field limit gives [ 16 ]: 

I(r )exp[-i/c(si - ^2) ■ r ]d 3 r' 
f I (r')d 3 r' 
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The beam splitter 22 was used for purposes of convenience 
only, and it is not a requirement of the instant invention that 
the detectors 10 and 12 be so oriented. The net effect, 
however, is to obtain signals from the detectors that are 
5 comparable to those that would have been obtained if the 
detectors had both been directly aimed at the sample 18 . 
Finally, in the text that follows the detectors will be spoken 
of as though they are both pointed directly toward the 
sample 18 , but those skilled in the art will understand that 
10 mirrors, beam splitters, etc. make it possible to position the 
detectors 10 and 12 in almost any arbitrary orientation. 

The input laser 14 is preferably positioned in the xy -plane 
at an angle $ (the “traditional” scattering angle) with respect 
to the x-axis. For 4>=90 degrees, the incident beam is thus 
15 parallel to the y-axis, and the x-axis preferably bisects the 
angle c|>, i.e., theta is the angular separation of the unit 

vectors s 1 and s 2 , which unit vectors lie in the x-z plane. 

Because (|), which is usually about 1 radian, is typically 
much larger than 0, which is preferably measured in milli- 
20 radians, the magnitude of the scattering wave vector is well 
approximated by the standard expression: 


In the previous equation, k=2jtn/X is the wave vector of 
the elastically scattered radiation having wavelength X in the 

medium with refractive index n. The unit vectors s ± and 

s 2 point from an origin in the scattering volume to the two 
detector positions. Note that the previous integration has 
been extended from a two dimensional surface integral to a 
three dimensional volume integral to accommodate scat- 
tered light coming from a three dimensional volume. The 

light intensity I( r ')is determined by the overlap of the 
incident beam, the multiple scattering within the sample 
volume, and the field of view of the detectors. A phase factor 
has been ignored in the above expression because it makes 

no contribution to |y( r i,r 2 ,0)| 2 , the quantity measured in the 
preferred embodiment of the instant invention. 

FIG. 1 shows a preferred two-detector scattering geom- 
etry for use with the instant invention, wherein the scattering 
angle § has been selected to be about equal to ninety degrees 
and the incident beam lies in the xy plane. Note that the 
geometry illustrated in that figure — including the choice of 
a particular value for the parameter § — is just one of many 
that could be used and it is well within the ability of those 
skilled in the art to devise other arrangements. The two 
detectors 10 , 12 in FIG. 1 are preferably oriented in the 
xz-plane and each forms an (extremely) small (e.g., mrad) 
angle 0/2 with respect to the x-axis. That being said, it is not 
essential to the operation of the instant invention that the 
detectors be oriented in any particular plane nor set at any 
particular scattering angle <|). It is anticipated, though, that 
generally the detectors and sample detection volume will be 
aligned in a plane which contains a direction that is perpen- 
dicular to the direction of the incident beam. The single 
scattered speckle is elongated compared to the multiple 
scattered speckle in these directions. It is required, though, 
that the two detectors 10 and 12 be aimed toward the sample 
in such a way that they detect substantially the same portion 
of the scattered field and that the angle between their 
respective lines-of-sight be very small. 

Further note in FIG. 1 that, strictly speaking, the detectors 
10 and 12 are not both aimed directly at the sample, but 
rather a beam splitter 22 allows detector 10 to be aimed in 
a direction that is roughly transverse to the direction in 
which detector 12 is aimed and still detect the scattered field. 


q=A%n sin(<j)/2)/A,, 

25 where “n” is the refractive index of the fluid suspension. The 
field of view of detector 10 is, thus, approximated by the 
following “gaussian tube”: 

I 1 (^>?,z)=exp(-a(y 2 +(z cos (0/2)-x sin (0/2)) 2 )). 

30 _ 

Here l/va is the radius of the 1/e point for the intensity if the 
field of view is uniformly illuminated. The field of view for 
detector 12 is similarly approximated by: 

^ I 2 (x;}?,z)=exp(-a(y 2 +(z cos (0/2) +x sin(0/2)) 2 )). 

The incident beam intensity is approximated by 

I ss (x,y, z)=Bexp(- p(z 2 +(x sin(<| >)+y cos(<|>)) 2 )), 

40 where the maximum intensity on-axis is given by the 
parameter “B”. The multiple scattered intensity is assumed 
to be uniformly distributed with magnitude “A” and it is 
further assumed to be incoherent with respect to the incident 
beam at the same location. This simply means that the 

45 incident and multiple scattered intensities are added 
together, rather than adding their field amplitudes and squar- 
ing the result. Finally, the finite size of the sample 18 is 
represented by a Gaussian cutoff function: 

50 I gc to,z)=exp(-25(x 2 +y 2 +z 2 )). 

Combining these factors together as they occur in the 
preferred measurement geometry yields: 

55 y ( fuh , o ) = 

exp(- &2kz sin(0/ 2)/i (x, y, z)h(x, y, z)I gc {x, y, z)(A + I ss (x, y, z))d 3 r 
f I l (x, y, z)/ 2 (x, y, z)I gc (x, y, Z)(A + I ss (x, y, z))d 3 r 

60 

Assuming a small focused laser beam and an intermediate 
detection width and a large sample volume, then the param- 
eters p, a, and 8 will be assumed to have the following 
relationship: p»a»8. Further, since the angle between the 

65 detectors, 0, is extremely small, the standard substitutions 
sin(0)«0 and cos(0)«l may be used. It is also fair to assume 
that (|)>>0 which leads, in turn, to the inequality: 
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(3 (sim|)) 2 »8. Given these simplifications and assumptions, 
a simplified expression for y(0) may be developed: 

^ exp {-k 2 6 2 j (8ar)) ^ exp (-k 2 6 2 / (4/5ar)) 

(® 2 (®fl 2 /4 + <5)) 1/2 ySV7si H0 5 

T( } “ A 2fl 

arV (5 £* sin0 

The square of this function gives the signal to noise ratio in 10 
equation 1 when t 1 =t 2 . The normalized intensity correlation 
function G (2 ) (t) is then related to the normalized field 
correlation function G (1 ) (t) by 

G< 2 \t)=1+y(0) 2 CG< 1 >(t)) 2 . 15 

Note that various assumptions and approximations have 
been made in the derivation of the foregoing that were 
designed to expedite the evaluation of some of the integral 
formulas. However, those skilled in the art know that these 

on 

substitutions were only representative of the many others 
that could have been made and are not critical to the practice 
of the instant invention. Thus, the various simplifying 
assumptions listed herein have been made for purposes of 
computational convenience only rather than out of any intent 
to limit the instant method to any one particular method of 25 
obtaining, solving, and evaluating the previous equations. 

PREFERRED APPARATUS 

In a preferred apparatus embodiment of the instant 3Q 
invention — and as is shown generally in FIG. 1 — a Helium 
Neon laser 14 (L=632.8 nm) is focused with a lens 16 of 
focal length 5 cm into a cylindrical sample cell 18 1 cm in 
diameter. The sample cell 18 is preferably immersed in a 
larger water bath cylinder 20 (6.5 cm diameter) to reduce 35 
parasitic scattering effects. Abeam waist radius of approxi- 
mately 30 /urn is calculated based on the beam size (0.68 mm 
radius for the intensity 1/e points) and Gaussian optics. This 
calculation ignores the influence of the index of refraction 
change of the water bath by using the focal length of 5 cm 4Q 
in air. The true beam diameter should be larger due to 
refraction effects. In practice, the instant inventor has found 
that a beam radius of 50 /um corresponding to p=4T0 s /m 2 
fits the experimental results best. 

As in conventional dynamic light scattering, the sample 45 
18 is illuminated by a tightly focused laser beam, such as is 
indicated at point 28 in FIG. 2, the scattering therefrom 
being captured by the two detectors 10 and 12 which are 
oriented at detector angle 0 with respect to each other. Note 
that, when the detector angle is zero, normal dynamic light 50 
scattering results and the recorded signal will typically be 
contaminated by a multiple scattered component. However, 
as the angular separation between the detectors increases the 
multiple scattering component of the signal decreases and 
ultimately only a single scattering signal remains. 55 

Two single mode fibers 40, 42 are preferably used as 
detectors and are aligned through a beam splitting cube 22 
to the focal spot of the laser inside the sample, beam splitter 
22 preferably being a polarization preserving beam splitter. 

A rough alignment is accomplished by reversing the detec- 60 
tion light path with a second auxiliary HeNe laser to arrive 
at an overlap of the incident and “detecting” beams. By way 
of example, in the preferred arrangement a HeNe beam 
directed through the fibers shows a divergence of 1.5 mrad. 

To determine the minimum detector angle necessary to 65 
record only the single scattering field, a scan over a range of 
detector angles 0 is performed. In the preferred embodiment, 
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the scattering angle, c|), is chosen to be 90° and the two 
detectors 10, 12 are arranged so as to deviate only a few 
mrads from (|) and are preferably positioned above and below 
the usual scattering plane. That being said, many other 
choices of scattering angle could have been used instead and 
that possibility has been specifically contemplated by the 
instant inventor. 

As is shown with more particularity in FIG. 2, the detector 
angle is preferably varied with a micrometer 26 that tilts the 
plane of one detector 10 with respect to the other 12. 
Although adjustment with a micrometer is the presently 
preferred method of changing the tilt angle and, hence, the 
angle between the two detectors 10 and 12, there are many 
other ways to accomplish this same task. By way of 
example, other ways that the detector angle might be 
adjusted include moving the two detectors laterally apart, or 
closer to (farther away from) the source, utilizing a curved 
beam splitter, etc. Thus, as used hereinafter, the term 
“micrometer” will be broadly construed to refer any mecha- 
nism that can be used to adjust the orientation of at least one 
of the detectors so as to change the detector angle 0 in very 
small (e.g., mrad) increments. 

According to another embodiment of the instant 
invention, an apparatus substantially as described above is 
used, but each detector 10 and 12 is further provided with a 
GRIN lens 41 that intercepts the speckle data before it enters 
the single mode fibers 40 and 42. As is well known to those 
skilled in the art, a GRIN (i.e., graded index) lens is a lens 
that has a refractive index that gradually changes from a 
relatively high value at its core to a lower value at its 
periphery. Thus, the refractive index decreases with the 
distance away from its core axis as a function of radius. In 
a GRIN lens, a process of continuous refraction bends rays 
of light such that SnelFs law is obeyed on a microscopic 
scale. Light ray containment now occurs by a process of 
continuous refraction. Among the many advantages of 
GRIN lenses are that they result in better detection effi- 
ciency. One GRIN lens that would be suitable for use with 
the instant invention is a model designed for 633 nm and 
sold by OZ-Optics, Canada. In the preferred embodiment, 
each GRIN lens 41 will be about 2 mm in length and made 
as an integral part of both ends of the single mode fibers 40 
and 42. 

If a GRIN lens is used, the detector field of view approxi- 
mated as a cylinder has a diameter of 1.1 mm (fiber 
divergence angle times fiber distance) and corresponds to 
a=3. 3 T0 6 /m 2 . A detector field cylinder diameter of 0.9 mm 
was found to improve the fit. By assuming that the multiple 
scattering sample volume is a sphere of «1 cm, the value of 
the parameter 8 can be taken to be 4T0 4 /m 2 . The cross 
correlation of the detector signals is preferably obtained with 
an ALV1-5000/E correlator (not shown). That being said, it 
is also possible to separately digitize the light signal 
received by each detector 10 and 12 and store those numeri- 
cal values for later computation. In that case, additional 
manipulations of the two signals would be possible 
including, by way of example, frequency domain analyses, 
time domain analyses, analyses of mathematical functions of 
the digitized detector signals, statistical analysis, cross cor- 
relations at other than zero lag, etc., any of which might 
potentially prove to be useful in a given situation. 

The tilt angle discussed previously that is created by the 
adjustments to the micrometer position 26 is preferably 
calibrated by reversing a laser beam through the single mode 
fiber 40 and noting the angular change on a screen placed 
some distance away from the apparatus. After the applica- 
tion of a given tilt, the detector typically has to be translated 
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slightly to align it with the incident beam and to maximize 
the detected intercept. It has been the inventor’s experience 
that this translation can be extremely sensitive to even a 
slight touch of the controls. This sensitivity is believed to be 
a principal reason for the scatter observed in the data points 5 
discussed hereinafter. 

According to still another aspect of the instant invention, 
there is provided an to apparatus substantially as described 
above, but wherein a polarizing filter is positioned to inter- 
cept the speckle data before that data reaches the detector. As 10 
is well known to those skilled in the art, a polarizing filter 
attenuates light passing therethrough that is not oscillating in 
a preferred plane. By placing a properly oriented polarizer 
between the sample and the detectors, some of the multiple 
scattering will be eliminated and the signal-to -noise ratio 15 
will be thereby increased. 

Finally, according to another useful embodiment of the 
instant invention, the apparatus was further arranged so as to 
be suitable for to measuring particle sizes in a polystyrene 
sample of 2% volume fraction. The arrangement of the 20 
detectors 10 and 12 is the same as that shown in FIG. 1 , but 
the sample 18 and water bath 20 are preferably replaced by 
a sample in a rectangular cell. The cell can be placed in an 
angular range of approximately 15-60° with respect to the 
incident laser direction (this angle is important for the 25 
scattering angle determination), and the laser is focused near 
the surface of the sample 18 . The detectors 10 and 12 are 
aligned to measure at the focal point and are preferably tilted 
about 1 mrad with respect to each other. This skin layer 
scattering technique works very well for the highly concen- 30 
trated volumes. It maximizes the visible single scattering 
while the tilted fibers suppress the attendant multiple scat- 
tering. 

ANALYSIS 35 

According to one preferred mode of data analysis, the 
intensity cross correlation function obtained from the two 
detectors 10 and 12 at a particular d etector angle 0' is fitted 
with a two-cumulant expansion of the form 40 

G^ 2 \T)=l+Y 2 (0')exp(-2MT+2vr 2 ) 

using the signal-to -background fitting factor y 2 (0’) defined 
previously, where T=t 1 -t 2 is the delay time, typically mea- 
sured in microseconds, and where G (2 ) (t) is the normalized 45 
intensity cross correlation function between the two detector 
signals at the fixed detector angle 0. Although this function 
might be estimated in many ways, in the preferred embodi- 
ment a non-linear least squares routine is used to produce 
estimates of the parameters u and v. Provided that the 50 
detector angle has been chosen in the range wherein single 
scattering predominates and wherein the multiply scattered 
signal is attenuated, the fitted value of the first cumulant u 
may then be used to estimate the particle radius r in the 
scattering medium through its connection with the Einstein 55 
diffusion coefficient D and the wave vector q: 
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Theoretically speaking, the previous equation requires that 
the particles whose size is to be estimated be freely diffusing 
and non-interacting as those terms are known in the art, e.g., 
particles in suspension. Of course, where there is a distri- 65 
bution of particle sizes — rather than a perfectly uniform 
sample — the instant method will return the z-average size of 


the particle distribution Thus, in the text that follows when 
a sample particle size is determined, it will be understood to 
be some sort of average or composite measure such as the 
z-average size if the sample contains particles of different 
sizes. 

In the previous equation, the parameter k B denotes the 
Boltzmann constant, T is the absolute temperature, n is the 
refractive index of the suspension, r\ the viscosity of the 
sample and L the laser wavelength. Additionally, the nor- 
malized second cumulant v/u 2 is an indicator of the “qual- 
ity” of the fit since, for absolutely monodisperse spherical 
particles in single scattering, the normalized second cumu- 
lant should vanish. The normalized second cumulant is also 
an indicator of the variance of particle sizes in the sample, 
i.e., for a sample containing particles with a distribution of 
sizes, this parameter can be an indicator of the polydispersity 
(i.e., the relative width of the size distribution). In the 
examples that are discussed below, normalized second 
cumulant values less then or equal to about 0.05 in the dilute 
single scattering limit are expected. 

EXAMPLES 

In the specific examples that are presented hereinafter, the 
sample consisted of polystyrene spheres of diameter 0.107 
jum manufactured by Duke Scientific (5.6% CV, Lot 16203) 
and diluted with water to 0.15 and 0.25% volume fraction. 
It should be noted that these examples are offered for 
purposes of illustration only and are not intended to limit the 
scope of application of the instant invention. 

FIG. 3 contains a plot of the signal-to -background for a 
range of different detector angles 0 for a lower volume 
fraction sample. The solid line that is drawn over the 
experimental measurement points is the expected behavior 
of a sample with an assumed ratio (A:B) of multiple to single 
scattering of 1:600. (PST 0.107 //m, 0=0.0015). Great care 
was taken in this particular experiment to measure y 2 near 
0=0. Note that for larger 0 the measured value of y 2 is more 
stable and reproducible and that the highest possible single 
scattering intensity auto -correlation function y 2 occurs at 
about 0.85. Therefore the fitting function (which is given by 
the square of the function for y(0) given previously) was 
multiplied by this factor. Additionally, note that measure- 
ments corresponding to negative values of 0 were also taken 
and that, as expected, those measurements proved to by 
symmetric about an angle of zero. 

FIG. 3 clearly indicates the presence of different speckle 
sizes: the central narrow peak corresponds to the multiple 
scattering speckle and the broader peak underneath is the 
single scattering speckle originating from the smaller single 
scattering volume. 

The previous tilt angle scan was repeated for an even 
more concentrated 0.107 /urn polystyrene sample, the results 
of which are presented in FIG. 4. There are two data sets in 
this figure, each of which was collected after passing the 
incoming speckle data through a polarizer: one data set was 
obtained with a polarizer that selected the same polarization 
(parallel) as the incident laser beam’s polarization plane, 
while the other set oriented to pass only the perpendicular 
component. The filled circles in FIG. 4 result from a two 
cumulant fit to the time dependent cross correlation function 
of the parallel polarization component, whereas the open 
circles represent the perpendicular polarization data. The 
solid lines on that plot have been calculated with an assumed 
ratio (A:B) of multiple to single scattering of 1:700 for the 
parallel and 1:0 for the perpendicular component. (PST 
0.107 /urn, 0=0.0025). 
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In FIG. 4, the parallel component shows a high intercept 
in the multiple scattering center region and then drops 
significantly to a shoulder region just as did the data in FIG. 

3. The intercept of the perpendicular polarization effectively 
vanishes for angles below-about -0.5 mrad, because the 5 
single scattering signal is severely attenuated by the polar- 
izing filter. The lines in this figure are calculated from the 
normalized intensity correlation function (G ( 2 ) (t), 
previously), but with an additional multiplying amplitude 
factor of 0.85 for the parallel and 0.65 for the crossed 10 
polarization signals as determined from the autocorrelation 
results discussed previously. A ratio of multiple to single 
scattering of 1:700 for the parallel and 1:0 for the perpen- 
dicular component is assumed: all other parameters 
remained unchanged. The single scattering contribution for 15 
the previous, more dilute, sample was smaller because the 
polarization of the incident laser was misaligned with 
respect to the scattering plane normal. Here, the incident 
polarization vector is perpendicular to the scattering plane. 

FIG. 5 contains a comparison of the results obtained 2 o 
during a test of a second 0.107 jum diameter polystyrene 
sample for the two polarization components. The filled 
circles are the radii obtained from a fit to the cross corre- 
lation function of the parallel polarization component, and 
the open circles are fits to the perpendicular polarization 2 s 
data. The cross hatched area is the expected size based on a 
single scattering analysis. (PST 0.107 pm, 0=0.0025). In 
FIG. 5, note that when the calculated fits are based on the 
perpendicular polarization component, the particle radii 
determined from the cumulant fit are substantially smaller 30 
than the actual value of 53.5 nm. This is primarily because, 
in this polarization orientation, only multiple scattering is 
observed. On the other hand, calculations based on the 
parallel polarization component gave very good particle size 
estimates where the angular separation was greater than 35 
about 0.5 mrad in magnitude. More particularly, at about 1 
mrad the influence of single scattering predominates and the 
calculations produce an acceptible size estimate. For detec- 
tor angle magnitudes less than 0.5 mrad, the multiple 
scattering component is not sufficiently attenuated and the 40 
particle size estimates are correspondingly affected. In 
general, as the angular tilt is increased the cross correlation 
loses most multiple scattering contributions and produces 
increasingly better estimates. 

FIGS. 3-5 confirm that the multiple scattering speckle 45 
decreases faster with detector separation angle then the 
single scattering speckle. It is always present regardless of 
choice of polarization orientation or the presence or absence 
of the single scattering speckle. The single scattering 
speckle, however, disappears when the detected polarization 50 
is perpendicular to the incident polarization (a forbidden 
“transition” for single scattering). The perpendicular corre- 
lation function decays faster than the parallel one because it 
contains relatively more multiple scattering. As the detector 
angle increases the correlation function begins to disappear. 55 
There is no detectable signal in the single scattering shoul- 
der. 

As was discussed previously, the value of the normalized 
second cumulant is a measure of the exponentiality of the 
correlation function. A small second cumulant (for mono- 60 
disperse particles) is a good indication of a single scattering 
correlation function. If multiple scattering adds a faster 
decaying contribution to the correlation function, the func- 
tional form will change. This non-exponentiality will sur- 
face in the form of a relatively large value of the normalized 65 
second cumulant. This contention is confirmed with the data 
in FIG. 6. As before, the filled circles are the fits to the cross 
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correlation function of the parallel polarization component, 
and the open circles represent the perpendicular polarization 
data. The cross hatched area represents the single scattering 
expectation for this sample. (PST 0.107 pm , <P=0.0025). In 
the multiple scattering region of the detector angle, the 
second cumulant is significantly above 10%. It drops to 
about 0.05 in the single scattering regime. 

CONCLUSIONS 

The present invention is thus an improvement over the 
well known ALV Dual Color light scattering apparatus in 
that it is significantly less expensive to set up and to 
maintain. Though the apparatus described above was origi- 
nally set up to examine the speckle structure as described by 
W. Meyer et al. in reference 13 of the following 
bibliography, its capabilities are enhanced relative prior 
techniques. Meyer et al. do not use GRIN lens fibers but 
utilize a sample temperature control bath and/or index 
matching vat as a lens for the fibers. In highly multiple 
scattering samples, where beam attenuation is significant, it 
will be best to detect the light scattering from a flat cell 
surface at the entrance point of the beam. This maximizes 
the single scattered signal and allows scattering angle stud- 
ies to be made. Under these practical situations, the GRIN 
lens system is superior. Furthermore, the GRIN lens system 
may be adjusted to optimize the single scattered signal by 
finding the angle 0 which is just outside the multiple 
scattering signal. This is not possible in the system of Meyer 
it al. 

The multiple scattering suppression techniques described 
herein will be useful in any dynamic light scattering experi- 
ment where multiple scattering is a problem and single 
scattering signals are desired, e.g., particle sizing in optically 
dense colloidal systems including paints, oil, gels, flocks, 
porous media, smoke, biomolecules; microsensors in moni- 
toring fluid quality and process control. 

While the invention has been described with a certain 
degree of particularity, it is understood that the invention is 
not limited to the embodiment(s) set for herein for purposes 
of exemplification, but is to be limited only by the scope of 
the attached claim or claims, including the full range of 
equivalency to which each element thereof is entitled. 
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What is claimed is: 

1. A method of determining a particle size using dynamic 40 
light scattering, wherein is provided a sample containing a 
plurality of particles, comprising the steps of: 

(a) directing a laser light beam at the sample, thereby 

creating a speckle field containing at least a multiply 
scattered signal and a singly scattered signal; 4 

(b) orienting a first detector so as to detect a predeter- 
mined portion of the speckle field; 

(c) orienting a second detector so as to detect substantially 
the same predetermined portion of the speckle field, 5Q 
said second detector being aligned at a predetermined 
small angle with respect to said first detector; 

(d) simultaneously detecting the predetermined portion of 
the speckle field within said first detector and said 
second detector, thereby creating a first detector signal 55 
and a corresponding second detector signal; 

(e) performing steps (b) to (d) at a plurality of different 

predetermined small angles, thereby producing a plu- 
rality of first detector signals and a corresponding 
plurality of second detector signals; and, 60 

(f) determining an estimate of the particle size using said 
plurality of first detector signals and said corresponding 
plurality of second detector signals. 

2 . A method according to claim 1 , wherein step (d) 

includes the step of: 65 

(dl) calculating a cross correlation between said first 
detector signal and said second detector signal, thereby 
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forming an intensity cross correlation value, and, 
wherein step (Q includes the step of: 

(fl) determining an estimate of the particle size using 
any intensity cross correlation values so calculated. 

3 . A method according to claim 1 , wherein step (f) 
includes the steps of 

(fl) selecting a particular first detector signal from among 
said plurality of first detector signals, 

(f2) selecting a corresponding second detector signal from 
among said corresponding plurality of second detector 
signals, 

(f3) cross correlating said selected first detector signal and 
said selected second detector signal, thereby producing 
an intensity cross correlation value, 

(f4) performing steps (fl) to (f3) for a plurality of different 
particular first detector signals and corresponding sec- 
ond detector signals, thereby producing a plurality of 
intensity cross correlation values, and, 

(f5) determining an estimate of the particle size using said 
plurality of intensity cross correlation values. 

4 . A method of determining a particle size using dynamic 
light scattering according to claim 3, wherein step (f5) 
includes the steps of: 

(1) determining from said plurality of intensity cross 
correlation values a particular detector angle, whereat 
the multiply scattered signal is attenuated with respect 
to the singly scattered signal, 

(2) orienting said first detector so as to detect a portion of 
the speckle field, 

(3) orienting said second detector so as to detect substan- 
tially a same portion of the speckle field, said second 
detector being aligned at said particular detector angle 
with respect to said first detector, 

(4) simultaneously detecting the speckle field within said 
first detector and said second detector, thereby creating 
a first detector signal and a second detector signal, 

(5) cross correlating said first detector signal and said 
second detector signal, thereby producing a cross cor- 
relation function, and, 

(6) determining an estimate of the particles size using at 
least a portion of said cross correlation function. 

5. A method of determining a particle size using dynamic 
light scattering according to claim 1, wherein step (f) 
includes the steps of: 

(fl) determining from said plurality of first detector 
signals and said corresponding plurality of second 
detector signals a particular detector angle, whereat the 
multiply scattered signal is attenuated with respect to 
the singly scattered signal, 

(f2) orienting said first detector so as to detect a portion 
of the speckle field, 

(f3) orienting said second detector so as to detect sub- 
stantially a same portion of the speckle field, said 
second detector being aligned at said particular detector 
angle with respect to said first detector, 

(f4) simultaneously detecting the speckle field within said 
first detector and said second detector, thereby creating 
a first detector signal and a second detector signal, 
(f5) cross correlating said first detector signal and said 
second detector signal, thereby producing a cross cor- 
relation function, and, 

(f6) determining an estimate of the particles size using at 
least said cross correlation function. 

6. A method of attenuating multiple scattering in dynamic 
light scattering measurements, wherein is provided a sample 
containing particles, comprising the steps of: 
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(a) directing a laser light beam at the sample, thereby 
creating a speckle field having a multiple scattering 
component and a single scattering component; 

(b) orienting a first detector to detect a predetermined 
portion of the speckle field; 

(c) orienting a second detector to detect substantially a 

same predetermined portion of the speckle field, said 
second detector being aligned at a predetermined small 
angle with respect to said first detector; 10 

(d) simultaneously detecting the predetermined portion of 
the speckle field within said first detector and said 
second detector, thereby creating a first detector signal 
and a second detector signal; 

(e) cross correlating said first detector signal and said 15 
second detector signal, thereby producing an intensity 
cross correlation value; 

(f) performing steps (b) to (e) at a plurality of different 
predetermined small angles, thereby producing a plu- 20 
rality of intensity cross correlation values; and, 

(g) determining from said plurality of intensity cross 
correlation values an optimizing detector angle 
between said first detector and said second detector, 
whereat said multiple scattering component is attenu- 25 
ated and whereat said single scattering component 
predominates. 

7. A method of attenuating multiple scattering in dynamic 
light scattering measurements according to claim 6, com- ^ 
prising the further steps of: 

(h) orienting said first detector to detect a particular 
portion of the speckle field; 

(i) orienting said second detector to detect substantially a 
same particular portion of the speckle field, said second 35 
detector being aligned at said optimizing detector angle 
with respect to said first detector; and, 

(j) simultaneously detecting the particular portion of the 
speckle field within said first detector and said second 4Q 
detector, thereby creating a first optimizing detector 
signal and a second optimizing detector signal at said 
optimizing detector angle. 

8. A method of attenuating multiple scattering in dynamic 
light scattering measurements according to claim 7, com- 45 
prising the further steps of: 

(k) determining an estimate of a size of said plurality of 
particles from said first optimizing detector signal and 
said second optimizing detector signal. 

9. A method of attenuating multiple scattering in dynamic 50 
light scattering measurements according to claim 8, wherein 
said particles are contained in a fluid suspension, wherein 
step (k) includes the steps of: 

(kl) cross correlating said first optimizing detector signal 55 
and said second optimizing detector signal, thereby 
producing an optimizing intensity cross correlation 
function, 

(k2) solving the following equation containing parameters 
u and v for at least the parameter u, 60 

G^ 2 \T)=l+Y 2 (0')exp(-2MT+2vr 2 ) 

where, y 2 (0') is a signal- to -background fitting factor at said 
optimizing detector angle, where t is a delay time, where 65 
G (2 ) (t) is said optimizing intensity cross correlation 
function, and where 0' is said optimizing detector angle, and, 


(k3) solving the following equation for the parameter r, 
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where r is an estimate of said size of said particles, where k B 
is Boltzmann’s constant, where T is an absolute temperature, 
where n is a refractive index of said fluid suspension, where 
r| is a viscosity of the sample, where X is a laser wavelength 
of said laser beam, and where c|) is said optimizing detector 
angle. 

10. A method of determining particle size in a sample, 
wherein is provided said optimizing detector angle of claim 
6 step (g), and wherein is provided a new sample, compris- 
ing: 

(a) directing a laser light beam at the new sample, thereby 
creating a new speckle field having a new multiple 
scattering component and a new single scattering com- 
ponent; 

(b) orienting a first detector to detect a predetermined 
portion of the new speckle field; 

(c) orienting a second detector to detect substantially a 
same predetermined portion of the new speckle field, 
said second detector being aligned at approximately 
said optimized detector angle with respect to said first 
detector; 

(d) simultaneously detecting the predetermined portion of 
the new speckle field within said first detector and said 
second detector, thereby creating a first detector signal 
and a second detector signal; 

(e) cross correlating said first detector signal and said 
second detector signal, thereby producing a new inten- 
sity cross correlation value; 

(f) determining a particle size of the new sample using at 
least said new intensity cross correlation value. 

11. A method of attenuating multiple scattering in 
dynamic light scattering measurements, wherein is provided 
a sample containing particles, comprising the steps of: 

(a) directing a laser light beam at the sample, thereby 
creating a speckle field having a multiple scattering 
component and a single scattering component; 

(b) orienting a first detector to detect a predetermined 
portion of the speckle field; 

(c) orienting a second detector to detect substantially a 
same predetermined portion of the speckle field, said 
second detector being aligned at a predetermined small 
angle with respect to said first detector; 

(d) simultaneously detecting the predetermined portion of 
the speckle field within said first detector and said 
second detector, thereby creating a first detector signal 
and a corresponding second detector signal; 

(e) performing steps (b) to (d) at a plurality of different 
predetermined small angles, thereby producing a plu- 
rality of first detector signals and a corresponding 
plurality of second detector signals; and, 

(h) determining from said plurality first detector signals 
and said corresponding plurality of second detector 
signals an optimizing detector angle between said first 
detector and said second detector, whereat said multiple 
scattering component is attenuated and whereat said 
single scattering component predominates. 

12. An apparatus for suppressing multiple scattering, 
wherein is provided a sample containing particles, compris- 
ing: 
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(a) a lasing apparatus, said lasing apparatus for producing 
a laser beam directed at the sample, thereby creating a 
speckle field signal when so directed; 

(b) a first detector positionable to detect at least a portion 
of the speckle field so created, said first detector 5 
generating a first electronic signal representative of any 
speckle field so detected; 

(c) a second detector positionable at a predetermined 
small angle with respect to said first detector and 
detecting substantially a same said at least a portion of 10 
the speckle field so created, said second detector gen- 
erating a second electronic signal representative of any 
speckle field so detected; and, 

(d) a mechanism for varying said predetermined angle 15 
between said first detector and said second detector in 
approximately mrad increments. 

13. An apparatus according to claim 12, wherein said 
mechanism for varying said predetermined angle between 
said first detector and said second detector is a micrometer. 

14. An apparatus according to claim 12, further compris- 
ing: 

(e) a beam splitter positionable between the sample and 
said first and second detectors, said first and second 
detectors detecting the speckle field through said beam 
splitter. 

15. An apparatus according to claim 12, further compris- 
ing: 
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(e) a polarizing lens between the sample and said first and 
second detectors, said first and second detectors detect- 
ing the speckle field through said polarizing lens. 

16. An apparatus according to claim 12, further compris- 
ing: 

(e) a correlator, said correlator being in electronic com- 
munication with said first detector and with said second 
detector, and said correlator cross correlating said first 
electronic signal with said second electronic signal, 
thereby producing at least a cross correlated intensity 
value. 

17. An apparatus according to claim 12 
wherein said first detector is a first single-mode optical 

fiber, said first single-mode optical fiber having a 
forward terminus positionable to detect at least a por- 
tion of the speckle field, and wherein said second 
detector is a second single-mode optical fiber, said 
second single -mode optical fiber being positionable to 
detect substantially a same said at least a portion of the 
speckle field. 

18. An apparatus according to claim 17 wherein said 
forward terminus of said first single-mode optical fiber 
further includes a GRIN lens and wherein said forward 

25 terminus of said second single-mode optical fiber further 
includes a GRIN lens. 



